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Abstract. A geometric formulation for stabilization of systems with one degree of undcractua- 
tion which fully solves the energy shaping problem for these system is given. The results show that 
any linearly controllable simple mechanical system with one degree of underactuation is stabilizable 
by energy shaping, possibly via a closed-loop metric which is not necessarily positive-definite. An 
example of a system with one degree of underactuation is provided for which the stabilization by 
energy shaping method is not achievable using a positive-definite closed-loop metric. 
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1. Introduction. One of the recent developments for stabilization of simple me- 
chanical systems is stabilization by energy shaping method. The central idea concerns 
the construction of a feedback for which the closed-loop system inherits the structure 
of a mechanical system. If such a feedback exists, the stability of the equilibria can 
be guaranteed if the Hessian of the closed-loop potential function is positive-definite. 
An important feature of the method, in case it is applicable, is providing a procedure 
that allows the design of nonlinear stabilizing feedbacks. 

The first classical appearance of the notion of potential energy shaping problem 
is in [34 . The investigation on the capabilities of this method continued in the 
Hamiltonian framework by looking at the properties of interconnected mechanical 
systems. The method is modified into the IDA-PBC method by introducing the 
notion of kinetic energy shaping [27]. An equivalent version of the IDA-PBC method 
in Lagrangian framework, so-called the Controlled Lagrangian method, is initiated 
by Bloch, Leonard, Marsden and Chang, [TO] [9] and its equivalence to the IDA-PBC 
has been proved in [14] [7j. In recent work, Chang, Woolsey and others have realized 
that the space of possible kinetic energy feedbacks can be enlarged by considering 
the addition of appropriate gyroscopic forcing (T2J [35] . It turns out that a necessary 
condition for stabilization of simple mechanical control systems by energy shaping 
method is linear controllability. For linear systems, linear controllability is also a 
sufficient condition for existence of a stabilizing feedback [36] [30]. In both methods, 
the question of energy shaping for a mechanical system reduces to solving a nonlinear 
system of partial differential equations. A large number of papers on energy shaping- 
method deals with finding a parametrization of solutions to this system of partial 
differential equations for a particular class of mechanical systems, for examples see [8] 
[37J EH [26]. 

A differential geometric approach to the kinetic energy shaping problem — the 
so-called A-method — has been presented in [SJ 2J [S] . A system of linear partial differ- 
ential equations is proposed for the kinetic energy shaping problem in terms of a new 
variable, A = G^G^, where G i and G c i are the open-loop and closed-loop metrics, 
respectively. The main idea of the A-method is that it transforms the set of quasi- 



* accepted for publication in SI AM Journal of Control and Optimization 

^Bahman Gharesifard is with the Department of Mechanical and Aerospace Engineering, Univer- 
sity of California San Diego, bgharesifardOucsd.edu 



1 



2 



BAHMAN GHARESIFARD 



linear equations for kinetic energy shaping into a set of overdetermined linear partial 
differential equations [3]. This method is extended to systems with gyroscopic forces, 
see [TJ. 

Lewis |24] has introduced an afhne differential geometric approach to energy shap- 
ing. In recent work, sufficient conditions for the existence of potential energy shaping 
are derived assuming that kinetic energy shaping has been performed |25] . The results 
are based on the integrability theory for linear partial differential equations developed 
by Goldschmidt [TB] and Spencer [33J . Gharesifard, et al. [TS] initiated a more system- 
atic geometric exploration of the kinetic energy shaping partial differential equations. 
The authors provide a set of sufficient conditions for the kinetic energy shaping. More- 
over, the kinetic energy shaping results are coupled with the integrability results of 
potential energy shaping |25j in order to provide a general approach for the total 
energy shaping. The technicalities of the geometric analysis of partial differential 
equations in [T5] and [23] might make it hard to comprehend the importance of these 
sufficient conditions. The results of the current paper fully relies on the integrability 
analysis of energy shaping partial differential equations and clarifies the importance 
of such results. 

Numerous systems considered in the literature on energy shaping have one degree 
of underactuation. In [3J the authors partially show that, under some conditions, 
systems with one degree of underactuation can be stabilized using energy shaping 
feedback. The results rely on a transformation of the system of partial differential 
equations and do not give any geometric insight into the energy shaping partial dif- 
ferential equations. In 115] the authors show that if E Q i is a simple mechanical control 
system with one degree of underactuation, for each bundle automorphism which sat- 
isfies the A-equation [6] 0] , there exists a closed-loop metric and a closed-loop poten- 
tial function which satisfy the energy shaping system of partial differential equations. 
Chang in [13] proves this result independently using the Cauchy-Kowalevski theorem. 
Note that this does not guarantee that such solutions are stabilizing ones. In fact, as 
we will demonstrate in this paper, the claim of |13] that the existence of such solutions 
to the kinetic energy shaping guarantees the existence of a positive-definite closed- 
loop metric which gives rise to an stabilizing energy shaping feedback is not true. 
However, in the current paper, we show that a linearly controllable simple mechanical 
control system with one degree of underactuation can be stabilized using an energy 
shaping feedback, possibly via a closed-loop metric which is not positive-definite. 

The paper is organized as follows. In Section [2] we recall the affine geometric 
formulation of the energy shaping problem. We give a summary of the integrabil- 
ity results for the partial differential equations in potential energy shaping in Sec- 
tion [3J We, directly and without proof, use the results of [T3J and [25]; thus we do 
not review the main integrability theorem of Goldschmidt [16]. A reader interested 
in understanding the formal integrability of partial differential equations is referred 
to [13 [III [111 [111 [331 EIJ EQl Ell [32] - Furthermore, we recall the A-method for kinetic 
energy shaping problem. Section [4] contains the main contribution of this paper: we 
show that all linearly controllable simple mechanical control systems with one degree 
of underactuation are stabilizable using energy shaping method. We fully character- 
ize a set of solutions to the kinetic energy shaping problem which is large enough to 
guarantee the stabilization by energy shaping method. 

Notation. The differential geometric notions used in modeling of simple mechanical 
systems are assumed here, and the unfamiliar reader is referred to PU El HI E31 E2] 



STABILIZATION OF SYSTEMS WITH ONE DEGREE OF UNDERACTUATION 



3 



for more details. The identity map for a set S is denoted by ids an d the image of 
a map / : S —> W by lm(/). For a vector space V the set of (r, s) -tensors on V is 
denoted by T£(V). By S&V and A^V we denote, respectively, the set of symmetric 
and skew-symmetric (0, fc)-tensors on V. Let A be a (0, 2)-tensor on V. We define the 
flat map A b : V — >• V* by (A b (it); v) = A(u, »),u,tieV, The inverse of the flat map is 
denoted by A* : V* -» V in case A b is invertible. We also define a similar notation for 
a (0, 3)-tensor A on V by 

(A b (u), w) = A(w, u, u), u, w G V. 

For S C V and W C V* we denote 

ann(S) = {a e V* | a(v) = 0, VoeS}, 
coann(W) = {v G V a(v) = 0, VaeW}. 

We denote by (E, tt, Q) a fibered manifold tt : E — > Q. The vertical bundle of the 
fibered manifold tt is the subbundle of Ttt : TE — > TQ given by Vtt = ker(T7r). We 
denote by Jfc7r the bundle of k-jets [3T]. A local section of tt is a pair (U, £), where Z7 
is an open submanifold of Q and £ is a map £ : ?7 — > E such that tt o £ = idy. If (£, £/) 
is an analytic local section of tt, we denote its fc-jet by jk£,- We denote an element 
of Jfe7r by jk£(x), where x & U. For more information about geometric properties 
of jet bundles see |3 1 j . A partial differential equation is a fibered submanifold Rfe C 
Jfc7r. Goldschmidt theorem investigates the conditions under which one can construct 
formal solutions of a given partial differential equation by constructing their Taylor 
series order by order |16j . 

2. Statement of the energy shaping problem. We assume that the reader 
is familiar with the affme geometric setup for simple mechanical systems [llj . A 
forced simple mechanical system is a quadruple S = (Q, G, V, J-" e ), where Q is an n- 
dimensional manifold called the configuration manifold, G is a Riemannian metric on 
Q, V is a function on the configuration manifold called the potential function and 
J- c : TQ — > T*Q is a bundle map over idq called the external force. We denote by V 
the covariant derivative with respect to the associated Levi-Civita connection. The 
governing equations for a forced simple mechanical system are 

V« (t) V(t) = -G« o dV(j(t)) + G».Fe(Y(t)), 

where 7 : / — > Q is an analytic curve on Q. 

Similarly, a simple mechanical control system is a quintuple S = (Q, G, V, T e , W) 
where Q is an n-dimensional manifold called the configuration manifold, G is a Rie- 
mannian metric on Q, V is a function on the configuration manifold called the potential 
function, J- c : TQ — > T*Q is a bundle map over idQ called the external force and W 
is a subbundle of T*Q called the control subbundle [TT]. The governing equations for 
a simple mechanical control system are 

V« (i)7 '(«) = -G« o dV( 7 (t)) + G» WW) + G^( 7 'W), 

where 7 : / — >■ Q is a curve on Q and it : TQ — > W is the control force. 

Given an open-loop simple mechanical control system S Q ] = (Q, G i, V \, J- Q \, W i), 
we seek a control force such that the closed-loop system is a forced simple mechanical 
system S c i = (Q, G c i, V c \, Fci), possibly with some external force. The reason for 
seeking this as the closed-loop system is that the stability analysis of the equilibria 
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for mechanical systems is well understood Chapter 6]. In this paper we assume 
that the open-loop external force J- Q \ is zero. 

It is well-known that presence of gyroscopic forces enlarges the space of possible 
closed- loop metrics [HI [24] while it does not change the total energy of the closed- loop 
system. As we see in next sections, systems with one degree of underactuation can 
be stabilized using energy shaping feedback without gyroscopic forces and thus our 
statement of energy shaping problem in this paper does not involve gyroscopic forces. 
In following, we present the statement of the energy shaping problem in the absence 
of gyroscopic forces. 

Definition 2.1 (Energy shaping problem in the absence of gyroscopic forces). 
Let S i = (Q, G i, V \, J~o\, VV i) be an open-loop simple mechanical control system 
with J- Q i = 0. If there exists a bundle map u s h p : TQ — > W i (called control) with 
u shp = — "kin — M pot such that the closed-loop system is a forced simple mechanical 
system S c i = (Q, G c i, V c i, 0) and 

1. Gii ° (*)) = V^ t)7 '(t) - V^ t)7 '(t), 

2- Upat(7(*)) = &oi &UvMt)) - dVMt)), 
then the control u s h P is called an energy shaping feedback. 

Throughout this work, we assume that the equilibrium point go £ Q is a regular 
point for W i- Moreover, we assume that the control codistribution W Q \ is integrable. 
This assumption is common in the literature and many examples fall into this case. 
The conditions of Definition 12 . 1 1 contain as unknowns the closed- loop metric G c i- One 
can observe that these equations involve the first jet of the unknowns. One can 
construct concretely a set of first-order partial differential equations as necessary and 
sufficient conditions for the existence of an energy shaping feedback. Let WcT'Q 
be a given subbundle and define the associated G i-orthogonal projection map P G 
r"(T*Q ® TQ) by 

Ker(P) = G^W. 

Note that P completely prescribes W. We apply P to the equation from part [TJ 
of Definition 12.11 to arrive at the following equation: 

^(V^ t) 7'W-V^ t)7 '(t))=0. 

Assume Q is an n-dimensional manifold and W is an integrable codistribution of 
dimension n — m. In adapted local coordinates the kinetic energy shaping partial 
differential equation is given by 

P? (^ll (^cl. I j.k + ^cUfcJ - Gcl^j,/) - G^Goi^fc + GoUfcj - Gol^jy)) = 0, (2.1) 

where i, j,k,l,r G {1, • • • ,n}, a G {1, • • • , m} and we denote the first derivative of 
GcUj with respect to q k by G c i y ^. Similarly, let P : T*Q — > T*Q/W i be the canonical 
projection on to the quotient vector bundle. We have 

P(G ol o G^Vd( 7 (*)) - dV i(j(t))) = 0. 

In local coordinates we have 

^(G oUj G cl 3 ^u - Voi, fc ) = 0, (2.2) 

where i,j, k G {1, • • • ,n}, a G {1, • ■ ■ , to} and we denoted the first derivative of V c i 
with respect to q k by V c i k- F° r more details on the affine differential geometric setup 
of energy shaping problem see [24) . 
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3. Summary of some integrability results. In [25] the potential energy shap- 
ing partial differential equation has been shown to be formally integrable under a 
surjectivity condition. An important corollary of this is that the choice of G c i affects 
the set of solutions that one might get for potential energy shaping. A bad choice 
of G c ; might make it impossible to find any potential energy shaping feedback. In a 
recent paper [15j , the authors show that the system of partial differential equations for 
kinetic energy shaping is formally integrable under a surjectivity condition. Moreover, 
they investigate the obstruction for integrability of the total energy shaping partial 
differential equations. 

Since the integrability conditions of the potential energy shaping partial differen- 
tial equations is an integral part of Theorem 14.31 we review this result in this section, 
without mentioning the proofs. Furthermore, we briefly recall the A-method for pro- 
ceeding with the kinetic energy shaping. We refer an interested reader to [15] for more 
details on the integrability character of the partial differential equations in A-method. 

3.1. Potential energy shaping. In this section, we explore aspects of potential 
energy shaping. We recall the results for potential energy shaping after kinetic energy 
shaping from [25] ■ Denote the bundle automorphism G b ol o GL by A c i. Define a 
codistribution W c i = A^ 1 (W i) and assume that this codistribution is integrable. 
Let PS = (Q x R, 7r, Q) be the trivial vector bundle over Q, so that a section of 
7r corresponds to a potential function via the formula q h4 (q,V(q)). We define a 
T*Q-valued differential operator Dd(V) = dV which induces a vector bundle map 
$ pot : Ji7r ->■ T*Q such that D d (V)(q) = $ P ot(ji V(q)). We denote by 

tt Wc1 : T*Q -> T*Q/W cI 

the canonical projection. 

Definition 3.1. Let E Q i = (Q, G i, V^i, J- \, W i) be an open-loop simple mechan- 
ical control system. The submanifold R pot C Ji7r defined by 

Rpot = {P 6 Jl7T | 7TW cl ° $ P ot(p) = nw c i ° A^dVoi} 

is called the potential energy shaping submanifold. One can easily observe that the 
"equation" representation of R po t is given by Equation (|2.2[) . 

Let 71"! : Jiir — > Q be the canonical projection. Lewis [25] gives a set of sufficient 
conditions under which the potential shaping problem has a solution. The proof 
follows from the integrability theory of partial differential equations; in particular, 
the potential energy shaping partial differential equation has an involutive symbol, 
see [20l El [H] for definition of involutivity. We recall the definition of (G i-G c i)- 
potential energy shaping feedback from [25] . 

Definition 3.2. A section TofWis called a (G i-G c i)-potential energy shaping 
feedback if there exists a function V c \ on Q such that 

T{q) - AddVd - dV oh q eQ. 

The following theorem establishes sufficient conditions for construction of a Taylor 
series solution to the potential energy shaping partial differential equation order-by- 
order. 

Theorem 3.3. Let Y, \ = (Q,G i,K,i, J^bWoi) be an analytic open-loop simple 
mechanical control system. Let G c i be a closed-loop analytic metric. Let po € R po t 
and let qo = ni(po) . Assume that qo is a regular point for W i and that W c i = A c ^ 1 W i 
is integrable in a neighborhood of qo ■ Then the following statements are equivalent: 
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1. there exists a neighborhood U of qg and an analytic (G i-G c i) -potential energy 
shaping feedback JeT" ( W) defined on U which satisfies 

* P ot(po) = AcjdV(gb) - d Voi (go) + A^dl^ao), 

for a solution V to R po t; 

2. there exists a neighborhood U of qo such that (i(A~ 1 1 dVoi)(g) G l2(W c i|g), where 
we denote l2(W c i| g ) = l(W G i| g ) fl A2(T*Q) and the algebraic ideal l(W c i|q) of 
A(T*Q) is generated by elements of the form 7 A uj with 7 G W c i| g . 

The theorem gives a set of compatibility conditions for the existence of a (G i- 
G c i)-potential energy shaping feedback. Moreover, one can give a full description of 
the set of achievable potential energy shaping feedbacks. Let a c i — dV v Let us 
use a coordinate system (q 1 , . . . , q n ) on U a neighborhood of q such that 

W c i| 90 =span(dg" l+ V-. ,dq n ). 

In these local coordinates we write the one form a c i as a c i — ctjdq 3 and compatibility 
conditions become: 

3.2. Kinetic energy shaping (the A-method). In following, we recall the 
so-called A-method in the absence of gyroscopic forces. The idea is to transform 
the kinetic energy shaping partial differential equations to an overdetermined linear 
partial differential equation, so-called the A-equation [6j [T2l [15] . 

Theorem 3.4. Let E„i = (Q, G i, V i, J~oi, VV i) be an open-loop simple mechanical 
control system. Let P G r"(T*Q(g>TQ) be the G \-orthogonal projection as above. Let 
Gci G L W (S^T*Q). 7/G o1 = G^oA for A G L w (T*Q<g)TQ), the following two conditions 
are equivalent: 

1. P(V% d X - V% al X) = , VATgL^TQ); 
faj V| o1 (Goi A) {PX, PY) = 0, and 

(b) V^ x G cl (Z, Z) + 2G C1 (V|° 1 APX, Z) = 2G i(v|°'PX, Z), 
where X,Y,Z G r u (TQ). 
For a complete version of the theorem and the proof, in the presence of gyroscopic 
forces, see [T2HH]- The set of A-equations have been proved to be formally integrablc 
under a surjectivity condition [T5] . 

3.3. An important corollary for systems with one degree of underac- 
tuaion. For systems with one degree of underactuation the potential energy shaping 
partial differential equations is always formally integrable. The main idea of the proof 
is that Equation (13. ip vanishes for m = 1, for details of the proof see [TB] . 

Theorem 3.5. If E Q i * s a simple mechanical control system with one degree 
of underactuation, for each bundle automorphism that satisfies the X-equation, there 
exists a closed-loop metric and a closed-loop potential function that satisfy the energy 
shaping partial differential equations. 

In the rest of these paper, we focus on stabilization of the closed- loop system. Ba- 
sically, we seek a solution to the energy shaping partial differential equation for which 
the Hessian of the closed-loop potential function can be guaranteed to be positive- 
definite. 
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4. Stabilization of systems with one degree of underactuation. In this 
section, we wish to determine the stabilizing solutions to the energy shaping partial 
differential equations for systems with one degree of underactuation. Throughout this 
section, let Q be an n-dimensional analytic manifold and E G i = (Q, G, V^i, W i) be an 
open-loop simple mechanical control system with one degree of underactuation. We 
denote the Hessian of a potential function V at qo G Q by Hess(V)(qo) G S2~T* o Q. In 
particular, we denote the Hessian of the open-loop potential function and the closed- 
loop potential function at the equilibrium point qo by Hess(V \)(qo) and Hess(Vd)(<7o)> 
respectively. 

Since the compatibility conditions of Theorem l3.3l are always satisfied for systems 
with one degree of underactuation, Theorem l3.5[ one can study the prolongation of the 
potential energy shaping partial differential equations instead of the original partial 
differential equations. Let (q 1 , . . . ,q n ) be local coordinates in a neighborhood U of 
qo G Q such that W i = spanjdq 2 , . . . , dq n } and let P be the projection of T*Q onto 
spanjdg 1 }. 

If we prolong the potential energy shaping partial differential equation and eval- 
uate the result at the origin, noting that dV c \(qo) = 0, we have 

P(G\qo)Gl 1 (q )d 2 V cl (v)(q ) - d 2 V ol (v)(q )) = 0, 
where v £ T go Q, i.e., 

G» 1 (< Zo )Hess b (F cl )(qo) - G»(<zo)Hes S b (Ki)(<Zo) - ^(q )(u\ qo ), (4.1) 

where u : TQ — > W i- If the system is linearly controllable, then one can design a 
control such that G B (go)Hess b (V \)(qo) + G$(qo)(u\ qo ) is diagonalizable and positive- 
definite. It is important to note that this does not necessarily imply that there exist 
G c i and V c \ such that G' 1 (qo)Hess l '(Vi = i)(go) is positive-definite, since the kinetic energy 
shaping partial differential equation puts restrictions on the achievable closed-loop 
metrics. However, we will show that, for systems with one degree of underactuation, 
the space of solutions of the kinetic energy shaping partial differential equations is 
large enough so that G'[(go)Hess b (T4i)(go) can be made positive-definite. We do this 
in the following steps. 

1. We first identify a simple class of solutions to the A-equation using Proposi- 
tion mu 

2. We show that this class of solutions is large enough to ensure that Equa- 
tion (14.11) holds with G|!j(<7o)Hess b (V^i)(go) diagonalizable and positive-definite. 

Let U be a neighborhood of the equilibrium point g £ Q an d let (q 1 , . . . , q n ) be local 
coordinates on U. In order to find the class of solutions mentioned in 1, we need to 
make some observations about the kinetic energy shaping partial differential equations 
for systems with one degree of underactuation. For these systems, the A-equation in 
the adapted local coordinate is given by 

-^ j :(& ll \\)~2S s kl G sl \\ =0, (4.2) 

where for i, j, k 6 {1, . . . , n}, are the Levi-Civita connection coefficients associated 
to G and i,k, s G {1, . . . , n}. Suppose we are seeking solutions to the A-equation that 
in local coordinates look like X(q) = Xjdq 1 ® where A- 6l and q G U, i.e., A is 
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constant. Then one can write Equation (|4.2[) as follows: 

- ^Ga) A} + - 2^ lGi2 ) A? + • • • + - 2Si lGin ) A" = 0. 

. (4 ' 3) 

Because >S is the Levi-Civita connection for G, the first term vanishes, leaving X\ 
arbitrary. One can rewrite Equation (|4.3p in the following fashion: 

n n 

^i^viA; (I. (4.4) 

i=l j=2 

where fe € {1, . . . ,n}. Thus, if X 3 2 — for j G {2, . . . ,n}, X(q) is a solution to the 
A-equation. Note that we further require that X(q) o G"(g) is symmetric. In the 
following, we describe the space of such solutions of the A-equation in an algebraic 
fashion. 

Let V be an n-dimensional K-vector space and let G £ S2V be a nondegenerate 
symmetric tensor. Let <E>g : V* ® V — >• A2V be the map defined by 

$g(A)(ui,i>2) = AoG(vi,v 2 ) -AoG(«2,«i), 

where i?i,«2 G V. The space of all tensors, A € V* (g> V, such that A o G is symmetric 
belongs to the kernel of <3?g and thus is of dimension " ( '" 2 +1 - ) , we denote this subspace 
by S G . 

Let {ei}" =1 be a basis for V and let {e l }™ =1 be its dual. Let W c V* be the vector 
subspace generated by {e 2 , . . . , e™} and denote its complement by E. We denote by 
S the space of all A € V* £g> V such that, if v £ coann(W), then A(v) £ coann(W), for 
all v £ V. A tensor A £ S can be written as 

n n 

A = Ale 1 (g> ei + ^ A^e 1 ® e 3 -, 

»=2 J=l 

where Aj S K and A] Gl for i £ {2, . . . , n} and j £ {1, . . . , n}. Thus the dimension 
of S is n(n — 1) + 1. If we denote the restriction of the map $g to S by $g |g : S — > A2V, 
then ker($cls) i s °f dimension n \ n ~ x > -\-\. If we additionally require that A £ kernel;;) 
be nondegenerate, we obtain a -dimensional subspace of V* g) V. 

Let Q be an n-dimensional analytic manifold and S Q i — (Qi G, V Q \, W i) be an 
open-loop simple mechanical control system with one degree of underactuation. Let 
U be a neighborhood of the equilibrium point q £ Q and let (q 1 , . . . ,q n ) be local 
coordinates on U such that W i| g = spanjeig 2 , . . . , dq n }, where q £ U. In following, 
we define a subspace of T*Q <E> T q Q which is large enough for stabilization of systems 
with one degree of underactuation. Consider the space of solutions to the A-equation 
that in local coordinates look like X(q) = X\dq l ® £ T*Q ® T g Q, where A^ G K 
and q £ U, and satisfies the followings 

1. X(q) o G'(g) is symmetric and nondegenerate; 

2. if v £ coann(span{ii(7 1 }) then X(v) £ coann(span{ii(7 1 }) for all v £ T q Q. 

We denote this subspace by 5? . The following proposition is a corollary of the algebraic 
discussion above. 

Proposition 4.1. 5? is an "^" 2 ~ 1 ^ -dimensional subspace of T*Q ® T g Q. 
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We wish to show that the space of solutions of the A-equation, described in Propo- 
sition 23 are large enough to guarantee that Gj 1 (<7o)Hess l '(V r c i)(go) can be made di- 
agonalizable and with positive real eigenvalues. If X(q) G then Equation (|4.1j) 
gives 

Hess\V cl )( qo )(JL, A) = JjW^feb)^, |j), (4.5) 
G,l x {q Q ){dq\dqi) = AiG»(<fo)(d?W), (4.6) 

where j 6 Asa result, we have the following proposition. 

Proposition 4.2. Let Q be an n- dimensional analytic manifold and S Q i = 
(Q, G, V i, VVoi) be an open-loop simple mechanical control system with one degree 
of underactuation. Let U be a neighborhood of the equilibrium point go G Q and let 
(q , . . . , q n ) be local coordinates on U such that W i| g = spanjdg 2 , . . . , dq n }, where 
q G U . Suppose that 

A = G»(go)rW(Ki)(?o) + G»(g )K ) 

is diagonalizable with real eigenvalues, where u\ qo : T 9o Q — > Woi|ij . Then there exists 
a closed-loop metric G c \ and a potential function V c \ such that 

1. G b = G\. y o A, where Ae^, 

2. G^ 1 (q , o)Hess b (V r c i)(qo) = A. 

Proof. We only need to show that if 1 holds, then G c i and V c \ can be selected so 
that 2 holds. Using Equations (|4.5I) and (|4.6p . we can write G c] [qo) in coordinates as 

/ \\a A}B A 
\\\& C J' 

where a G R, B G L(M"" 1 ,M), and C G S 2 R"~ 1 are such that a = G^(dq 1 ,dq 1 ) and 
B(dq\dq j ) = G^(dq 1 ,dq j ) for all j G {2,...,n}. Similarly, Hess b (V c i)(qo) can be 
written as 




wherefc G R, B G L(R"- 1 ,R), andC G SsR"" 1 are such that k = Hess b (Ki)(<to)(3§r, afr) 
and B( * » ) = Hess b (Ki)(go)(«fr, ^j) for all j G {2, . . . ,n}. Thus we have 



G» 1 (g )Hess b (Ve 1 )((7o) = G^Hess^Xgo) + 




where 

1. Li = kB T + jtCB t G L(R,R™- 1 ) and 

2. L 2 = A}B23 T + CC T G L(R ,l_1 x R"- 1 ) 

can be set to any value by appropriate choice of C and C. □ 

Theorem 4.3. Let £ i = (Q, G i, V i, W i) be a linearly controllable open-loop 
simple mechanical control system with one degree of underactuation and with qo G Q 
an equilibrium point. Then the system is stabilizable at qo using an energy shaping 
feedback. 
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Proof. The integrability of the energy shaping partial differential equations en- 
sures that formal solutions exist. Furthermore, Theorem 13.51 implies that the ob- 
structions of Theorem 13.31 are satisfied for systems with one degree of underactua- 
tion. If the system is linearly controllable, then one can design a control such that 
G"(<7o)Hess l '(T4i)(qo) + G"(go)(u|g ) * s diagonalizable and positive-definite. Proposi- 
tion then guarantees that G c i can be found such that it satisfies the kinetic energy 
shaping partial differential equations, by choosing A 6 S", and taking 

G» 1 ( (Zo )Hess b (K 1 )(go) = G«(g )Hess b (Ki)(go) + G»(g )« ) 

to be diagonalizable with positive real eigenvalues. □ 

Note that this proof does not require that the closed-loop metric be positive- 
definite and in fact, there are cases for which energy shaping is not possible with 
positive-definite closed-loop metrics; an example of this is presented in Example [SJ 
The following proposition clarifies when it is necessary to perform kinetic energy 
shaping for systems with one degree of underactuation. 

Proposition 4.4. Let Q be an n- dimensional manifold and let S Q i = (Q, G, V a \, W i) 
be a linearly controllable simple mechanical system. Let U be a neighborhood of qo 6 Q 
such that Woi = spanjdg 2 , . . . ,dq n }. If Hess(K,i)(^|r, g|r) > 0, the system can be sta- 
bilized around its equilibrium point qo without kinetic energy shaping. 

Proof. We shall show that S i is stabilizable using an energy shaping feedback 
with G c i = G. Equation (|4.1I) then reads 

Hess\V cl )(q Q ) = Hess b (V ol )(q ) + u\ qo , 

where u is a feedback. Note that since Hess(Vci) is symmetric, it is positive-definite 
if and only if all of its principal minors are positive. The first principal minor of 
Hess b (Vci) is positive. Then, by linear controllability, one can choose the controls so 
that the system is stabilizable at the equilibrium point go, similar to Proposition 14.21 
□ 

Next, we present an example of energy shaping for simple mechanical systems 
with one degree of underactuation for which the energy shaping is possible only via 
a closed-loop metric that is not positive-definite. 

5. Example. Consider the stabilization problem for a simple mechanical control 
system £ = (R 2 , G, V u 0, W ol ) at the origin q = € R 2 , where 

1. G = ((q 2 ) 2 + l)dq x ® dq 2 + {{q 1 ) 2 + l)dq 2 <g) dq 2 , 

2. Vol --(g 1 )^ 2q 1 q 2 + (q 2 ) 2 , and 

3. Woi = spanjdg 2 }. 

This system is linearly controllable at the origin. We show that, for any solution of 
the A-equation, the constant term in the Taylor expansion of X\ is always zero. In 
order to show this, we need to modify Equation (|4.4|) by adding an extra term, since 
A, in a neighborhood of go, is not necessarily chosen from S*. We have 

E( G « JT + E(^ G *i - ^G^Aj) = 0, 

i=l q j=2 
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for all k G {1, ...,n}. For this example, by substituting the nonzero Christoffel 
symbols, we have 



((9 2 ) 2 + l)^f + 2g 2 A? = 0, (5.1) 
((9 2 ) 2 + l)ft-VA? = 0. (5.2) 



It is clear that \\{qo) can be chosen arbitrarily. Consider formal expressions for \\ 
and A}: 

\\ = Coo + Cm 1 + Coiq 2 + C 20 (g 1 ) 2 + C 02 (q 2 ) 2 + C xl q\ 2 + • • • , 
\\ = D o + D 10 q 1 + D 01 q 2 + D 20 (q 1 ) 2 + D 02 (q 2 ) 2 + D 11 q 1 q 2 + • • • , 

where dj, Dij G K for i,j G Z>o- If A} and \\ satisfy Equations (|5.ip and (|5.2[) . 
then Cn = Doo = 0, i.e., A^(go) = 0- Thus the closed-loop metric at the origin has 
the form G c i(<Zo) = ^dq 1 ® dg 1 + ^rfg 2 <& dq 2 , where a,c G R\{0} and A}(g ) = a- 
Equation (|4.5[) implies that 



Hess b (^ cl )(go) = 



zl 1 

a a 

1 k 



where tel. Thus 



G» 1 (< Zo )Hess b (Ki)(9o) = ( 2c 2 l k 

\ a 

It is easy to see that one has to choose ^ < and ck > 2 in order to make 
Gci(<?o)Hess l, (V^i)(q , o) positive-definite, i.e., none of the achievable closed-loop metrics 
is positive-definite. However, one can choose a, c, k G R so that G^(go)Hess b (Ki)(go) 
is positive-definite, for example a = — c = || , and fc = 1. 
Remark. If we take the open-loop metric given by 

G = ((? 2 ) 2 + ljdg 1 8> d? 2 + ((q 1 ) 2 + l)dq 2 ® dq 2 + 2q 1 q 2 (dq 1 ® dq 2 + dq 2 ® dg 1 ), 

then A 2 (go) need not be zero and the system can be shown to be stabilizable by the 
energy shaping method with a positive-definite closed-loop metric. This reveals that 
a slight change in the structure of the open-loop Levi-Civita connection has a huge 
impact on the achievable closed-loop metrics. 

6. Summary. In this paper, we fully solved the problem of stabilization of sys- 
tems with one degree of underactuation. The result completely relies on the inte- 
grability analysis of partial differential equations involved in energy shaping. We 
illustrated that all linearly controllable simple mechanical control systems with one 
degree of underactuation can be stabilized using an energy shaping feedback, with 
closed-loop metrics which are not necessarily positive-definite. We also characterized 
the simple mechanical systems for which the energy shaping is achievable without 
kinetic energy shaping. Finally, we gave an example of a simple mechanical control 
system with one degree of underactuation for which there exists no solution to the 
energy shaping problem with positive-definite closed-loop metric. The results give 
some useful insight about the structure of kinetic energy shaping. 
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